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1 Introduction 



Recently, the study of vr-junctions with an intrinsic phase difference = vr in the ground state 
has gained much interest. Many experiments on high-Tc superconducting weak hnks have been 
interpreted in terms of vr-junctions proposed for the superconductors with d-wave pairing 
For conventional superconductors with s— wave pairing, vr-junctions have been suggested by 
Bulaevskii et al. [0] to arise in presence of magnetic impurities in the barrier, in connection 
with spin flip assisted coherent tunneling. Also, vr— coupling was predicted theoretically to 
exist in superconductor-ferromagnet {S/F) multilayers by Radovic et al. 0. The evidence for 
vr- coupling, manifested via characteristic oscillations of the superconducting critical temper- 
ature Tc with magnetic layer thickness, was sought experimentally in several superconductor- 
ferromagnet 0, 1^, 1^ and superconductor-spin glass systems . 

Other manifestations of nontrivial coupling ^ 0) could be found in the density of states 
(DOS) and, in particular, in the appearence of the Andreev bound states at the Fermi level. 
Similarly as in high-Tc superconductors, these zero-energy states (ZES) should be manifested 
as the zero-bias conductance peaks (ZBCP) in tunneling spectroscopy measurements |p. Re- 
cently, Andreev reflection has been studied theoretically in S/ F tunnel junctions by de Jong and 
Beenakker |^, in 5'/F/S' junctions by Kuplevakhskii and Fal'ko |]TU| and by Tanaka and Kashi- 
waya [|TT|, and in atomic-scale S/F superlattices by Prokic and al. [|1^. Recent experiments 



T3| , |n| have shown that Andreev reflection can be strongly suppressed by spin polarization in 



the ferromagnet. This is the case for point contacts with transparent barriers, whereas the spin 
polarization can enhance the Andreev reflection in the S/F/S junctions with sizeable interface 



scattering 15 



In the previous papers we calculated DOS and spontaneous currents in Josephson S/F/S point 
contacts [|^, [l^, considering the combined effect of d— wave pairing in S electrodes and of the 



exchange field in F barrier. 

In the present paper we calculate the quasiparticle density of states in conventional supercon- 
ductor Josephson weak links with a transparent ferromagnetic metal barrier. The results are 
simpler than in the general anisotropic case |T^, due to the isotropy of the order parameter in 
S. The phase difference at the link, related to the flow of the supercurrent, in presence of the 
exchange energy h in F acquires an additional, "magnetic" contribution. We study the change 
in the quasiparticle energy spectrum induced by and h and discuss the appearence of ZES 
and SPGS. 

The paper is organized as follows: in Section II we present briefly the results of the quasiclassical 
theory of superconductivity for S/F/S weak links: the quasiparticles Green's functions in the 
barrier, the resulting densities of states and the Josephson supercurrent. We show how the local 
densities of states can be expressed in terms of energies of bound states and give the equations 
for the bound states, from which follow the conditions for ZES and SPGS. We calculate the 
supercurrent through the link and discuss its relationship with Andreev bound states. Section 
III contains a discussion of the numerical results, including the comparison with the case of the 
normal metal (A^) barrier, h = 0, and a brief conclusion. 



2 Quasiclassical Theory 

An efficient method for calculating local spectral properties of superconductors is the quasi- 
classical theory of superconductivity |TB[. In the clean limit, and assuming that the magnetic 
influence on superconductivity is limited to that of the exchange energy in the barrier, one can 



use the Eilenberger quasiclassical equations in the presence of the exchange field [|19|, ^ 
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2.1 Quasiparticles Green's function 



In the previous papers we have solved the quasiclassical equations for an S/F/S weak hnk, with 
a transparent, thin and short F barrier of thickness 2d, consisting of a monodomain ferromag- 
netic metal with constant exchange energy h. Here we present the results for the conventional 
case of s— wave pairing in S electrodes. Both 5* and F metals are assumed clean, with same 
dispersion relations and with same Fermi velocity vq (electron scattering on impurities in S can 
be neglected if I ^ C,o, where / is the electron mean free path and superconducting coherence 
length, and if /i ^ hvo/l in F). 

Taking the x-axis perpendicular to the barrier (Fig. 1), we assume a step-function variation of 
the pair potential 



A{x) =AQ{-d-x)+AQ{x-d). (1) 

where A = A(T). In principle, the pair potential should be determined self-consistently. 
However, this would greatly complicate the calculations without shedding much light on the 



exchange- interaction- induced states in the S/F/S case [|10| , [2l[] , which is the main subject of 
this paper. In the S/N/S case, the validity of our model requires a weak proximity effect in S, 
as it may be the case at low T and for thin and short barriers. 
The quasiclassical Green's function in the barrier is 

cj„cos7„ + zfi„sin7„ 

9 = 9l = 7^ —■ : (2) 

iln cos7„ + lujn sm 7„ 



where Un = iTT{2n + 1) are the Matsubara frequencies {h = ks = 1), = \l^n + I^P^ 

2hd 2iuJnd , . 

7n = 7T + 3 

/ Vq cos if Vq COS (f 

and if is the angle between the direction vq and the x— axis |TB|, |T^. 
Outside the barrier, > d, we find 

where a„ = 2VLn/vQ cos if and a =1, J, denotes the spin orientations with respect to the direction 
of the exchange field. 

For the opposite spin direction, the corresponding Green's function is obtained by changing 
h^-h 

g^(h)=g^i-h). (5) 

2.2 Densities of states 

The quasiparticle spectrum follows from the retarded Green's functions, obtained by the ana- 
lytical continuation of gi{h) and gi{—h). With our definition of g^, the partial density of states 
(PDOS), which is the angle-resolved DOS, is given by 

N„{x, ip, E) = lim3fJ5f^(a;, (p, iun E + iS), (6) 

0— >0 

where E is the quasiparticle energy measured from the Fermi level. 
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DOS is obtained by averaging PDOS over the angle assuming spherical Fermi surface 

N„{x,E)/No = r^\<fsm<fV{if)N,{x,if,E), (7) 



where A^o = Tnkp/2T['^, and T>{ip) is the normalized barrier transmission probability |23|. We 
model the barrier with a uniform probability distribution "D^ip) = 1/ simpdip, within an 
acceptance cone of angle 2^9^ about the interface normal, and zero outside the cone. The 
relation between PDOS in S: 

N^ix, if, E) = \imlkgl{x, if, luJn ^ E + i5), (8) 

and in the barrier, Eq. (||), is obtained using Eq. (^). For \E\ < \A\ Eq. (||) reduces to the 
simple relation 

iVf(x,<^,E) = e°('^±^)iV,(<^,E) , (9) 



where a = 2Vt/vQ cos(^, and Vt = y |Ap — E"^. The influence of the barrier is seen only in the 
close vicinity of the interface. Deep in S we obtain, as expected, the well known bulk result 

N'{E)J-^Q{\E\-\/\\). (10) 

Since PDOS in the electrodes and in the barrier are related, we consider in the following only 
the position-independent PDOS and DOS in the barrier. 

Measuring all energies in units of A, E = E j IS., = fi/A, h = hj IS, for PDOS we get 

iVi(<^,|E|<l) = i5(E + ficot7), (11) 
where 8 is the Dirac delta function, and 

iV,(y,|^|>l)= f '^'^^ , (12) 
E'^ — cos"' 7 

where 7 = </)/2 — d{E — h) / cosip, d = 2dA/vo = {2/n)d/^o, and ^0 = vq/ttA. 
Bound states exist only for energies l-El < 1 and satisfy the equation 

tan7 = -£ , (13) 

corresponding to the poles in the analytical continuation of the Green's function, Eq.(|^). A 
more convenient form of Eq. ([13|) is 

sin 7 = ±Q, cos 7 = =F-E (14) 

with the condition 

sign(_E) sin7 C0S7 < . (15) 

Note that the bound states with opposite spin orientation, spin up, are given by Eq. (|T^) (or 
by Eqs. (p!4D, (|15D) with h — /i. There is a spin splitting of bound states. 
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Another consequence of Eq. (|T^) is that bound states at the Fermi level (ZES) appear for 
cos 7 = 0. For = this gives 



2 ' '2 

Calculating DOS for = 'n /2, we find 



-±dh= {2k + l)-, A; = 0,±1,±2,... . (16) 



< l)/iVo = (17) 
2Trd\E-h\ £ ^^"^^ 



~\ 2 



(^0/2 + nvr + arctan Vt/ E 

where A = sign{E — h) ($ + mr) and $ = aictanQ/ E + (j)/2 — d{E — h) . Let us denote by Ei the 
zeros of Eq. (pTsP for (f = 0, which can be written in the form $(-Ej) + njvr = 0, nj = 0, ±1, .... 
Then, Eq. (|17|) can be rewritten (see Appendix A) in the form 



Nm < l)/iV„ = 2.dlE -M{t^- - + - ^'^'^ - 



i=l 



(^0/2 + n jTT + arctan 



~\ 2 



^ e(^-E) - e(^ - h) 

n=-oo (0/2 + nvr + arctan n=n+ (0/2 + nvr + arctanil/E' 



where m is the total number of bound states E^, and integers ra" are defined by the 
conditions $ + nvr > for n > and $ + nvr < for n < n~, respectively. Integers 
rii G [n^ + 1, n"*" — 1] are chosen so that sign ($(£') + njTT^ = sign (^E — hj. 
The physical meaning of this result is the following: due to the isotropy of the order parameter 
the averaging over the quasiparticle propagation angle (p results in an average over one-chanel 
weak links of different lengths. Therefore there is no qualitative difference between DOS and 
PDOS, the peaks in DOS corresponding to peaks in PDOS for 93 = situated at the energies 
E,. 

For normal metal barrier, h = 0, our result agrees with Refs. [24] and [22] with small corrections 
of DOS for E > Ei. For = the total number m of bound states depends on the thickness. 



m = 2{[d/'K] + 1), where the symbol [A] represents the greatest integer less than A [|22|. In this 
case, ZES may appear only when the phase difference at the link is different from zero, Eq. 

For ferromagnetic metal barrier, h ^ 0, the ground state phase difference (pgg has been calculated 
T7| by minimizing the energy of the link, and it was found that (pgg = for weak barrier 
influence, 2hd < 1, and (pgs = vr for stronger influence, 1 < 2hd < 4. Larger values of 2hd would 



correspond to the decoupling of S electrodes [0. Combining this result with the condition 
for ZES, Eq. (p!6D, one can easily see that ZES can occur for 7^ (f)gs only. Due to the spin 
splitting, SPGS may appear [|l^, |ri|, Analyzing Eq. (|T3D for = 0, we find that for d > 7r/2 
there is always at least one bound state for each spin orientation in the interval —1<E<0. 
Therefore, the spin polarized ground state, where only one spin orientation negative energy 
levels are populated, may exist only in thin barriers, d < tt/2. They actually appear (see 
Appendix B), with spin up orientation, provided that (f)/2 + dh belongs to the (open) interval 
({2k + l)vr/2, {k + l)7r — dj, whereas 0/2 — dh lies outside the above interval, and vice versa 
for spin down orientation. 

We emphasize that all above analytical results are in complete agreement with numerical cal- 
culations of PDOS and DOS, which we performed directly from Eqs. (0) and (P). In numerical 
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calculations, we have everywhere taken = 7r/2. For smaller there is no qualitative differ- 
ence in DOS , the number and the positions of the peaks are the same, the relevant direction 
being that perpendicular to the barrier, = 0. This is not the case for c?— wave pairing in 
cuprates, as shown in Ref. [llB]. 



2.3 Supercurrents 



The supercurrent trough the S/F/S* junction is also obtained from the Green's functions in the 
barrier, Eqs. (||) and 



j = -2^e7riVoT^(vo 



where (■ ■ ■) denotes the angular averaging over the Fermi surface. Putting u = 1/ cosip we get 



4A 



EE 



du, cr =1, 1 . 



(19) 



This gives the supercurrent through the barrier of the area S, I = jS as a function of (p, 
temperature T, and of parameters measuring the influence of the F barrier, h and d. The 
normalizing current is Jq = 2A/eRN, and the normal resistance is given by R]^^ = c^voNqS, 
where Nq is the density of states at the Fermi surface. 

From Eq. (|19D we see that poles of ga determine both the quasiparticle energy spectrum and 
the phase and exchange field dependence of the supercurrent. That the supercurrent is carried 
by the Andreev bound state is particularly clearly seen at T = 0. Taking T ~^ I duj/271, it 
is easy to show that Eq. (|19D at zero temperature can be rewritten in the form 



/(0) = ^/oE 



du 



1 u-" 



(20) 



where Ej^i and Ej^^ are negative energy levels, obtained from Eq. (|l^) with 7 = 7(±/i). Since 
the pair potential is assumed to be spatially constant in the S electrodes, the free energy of the 
junction is related to the supercurrent by the well-known relation 



/(0) = 2e 



OF 



(21) 



and 



2eF(0) = -/oE / ^[E,,^ + E,,^]. 



1 



[22) 



In particular, in the limit of strong exchange field, ^ 1, we obtain 



IT du 



sin(^ + dhu)sign ^cos(^ -|- dhu)^ + sin(^ — dhu)sign ^cos(^ — dhu)^ 



and 



TT r°° du 
2ef (« = --/„/ 



cos(— + dhu) 



+ 



cos(— — dhu) 



(23) 



(24) 



Therefore, even in this limit, the supercurrent is carried by the Andreev bound states of both 
spin orientations, and the phase dependence is not sinusoidal, in contrast to the results of 



Tanaka and Kashiwaya for ferromagnetic tunnel junctions [jTT|. 
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3 Summary and discussion 



In this Section, we discuss the results obtained for S/N/S and S/F/S weak hnks, with thin 
metalhc barrier, 2(i/^o < (c? < 1). 

The results of numerical calculation of the supercurrent at T = as a function of (j), for some 
characteristic values of h, are shown in Fig. 2, and the results for PDOS and DOS in Figs. 3-6. 
In all examples we have taken 2d = C,o {d = I/tt). 

In both cases of normal {h = 0) or ferromagnetic {h ^ 0) metal barrier, the quasiparticle 
spectrum in the barrier is gapless, with characteristic peaks inside the superconducting gap, 
\E\ < A. Due to the isotropy of the order parameter in S, there is no qualitative difference 
between PDOS for ^9 = and DOS. The peaks in PDOS, corresponding to the bound states, 
are reflected in DOS as peaks placed at the same energies. This is illustrated for /i = in Fig. 
3, but holds for /i 7^ as well. 

For thin normal metal barrier, the number of bound states is m = 2 or m = 1. The spectra are 
strongly influenced by the phase difference 0. A single peak, corresponding to ZES, appears for 
= TT at the maximum of the supercurrent through the link. Fig. 2(a). In the ground state, 
= (pgs = 0, DOS starts from zero at ii^ = and has two peaks at E ^ 0, Fig. 3. 
For ferromagnetic metal barrier, h 0, the bound states for two spin orientations (spin up and 
spin down with respect to the exchange field orientation) are no more degenerate. However, 
both kinds of bound states contribute to the supercurrent through the link. The dependence 
/(0) is shown in Fig. 2(b) for (pgs = and for (pgg = vr. DOS curves are splitted with respect 
to the case h = 0, and the number of peaks is doubled. An example with = is shown in 
Fig. 4. 

Now, the conditions for ZES, Eq. ([T6|) , are ± 2dh = vr. The main effect of /i 7^ is to induce 
an additional, "magnetic" phase difference — » ± 2dh, for the electrons traveling in the 
direction (f = 0. Thus, the formation of ZES is due to an appropriate value of the total phase 
difference. When a ZES appears for a given value of and h for one spin direction, there is no 
ZES for the same for the opposite spin direction, h —h, except if = and dh = tt/2. 
The latter case is illustrated in Fig. 5 for = 7r^/2. Here (pgs = tt, and ZES is obtained outside 
the ground state, for = 0. 

Appearence of spin polarized states, with a single spin orientation at negative energies, is 
illustrated in Fig. 6, for h = 1.5. For given h and d, the spin orientation of the polarized state 
strongly depends on the phase difference 0. In the given example, the negative energy state 
in the ground state, = (pgs = 0, has the spin down orientation, and for = tt, the spin up 
orientation, in accordance with the condition ( p7|) in the Appendix B. 

In conclusion, we have analyzed the quasiparticle spectra in clean S/ F/ S contacts, with s— wave 
pairing in S electrodes. In this case DOS for each spin orientation is asymmetrical, and strongly 
exchange field and phase dependent. For both ferromagnetic and normal metal barriers we have 
found conditions for the formation of zero energy bound states, and in the case of ferromagnetic 
barrier, for the formation of spin polarized ground states. In all cases, zero energy bound states 
appear only outside the weak link ground state. 

Appendix A 

In this Appendix we show how in the case of s— wave pairing, DOS for E < 1, Eq. ([T^), can 
be expressed in terms of bound states energies, Eq. (^). Since the function $ = arctanr2/i^ + 
0/2 — d{E — h) is bounded in the interval \E\ < 1, there is a minimum integer n"*" such that 
$min + n'^T^ ^ 0. This mcaus that for each n > and each i?, $ + titt > 0. Similarly, there 
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is a maximum integer n such that for each n < n , and for each $ + nvr < 0. Thus, Eq. 
can be written in the form 

N,{\E\<l)/No = 2nd\E-h\ \ - E) 

I n=-oo ( 0/2 + nTT + arctan n/E) 

, t , , t 7 , J , (25) 

n=n+ I (f)/2 + mr + aictanQ/ Ej n=n-+i I (f)/2 + mr + aictan Q/Ej 

where A = sign(£' — /i) ($ + nvr). In this expression the first two sums are continuous func- 
tions in the whole interval \E\ < 1. In the third sum, all terms for which the condition 
sign (^^{E) + nvrj = sign (^E — hj is not satisfied are missing, giving rise to finite discontinu- 
ities. These discontinuities occur at bound states energies E = Ei ior ip = 0, as can be seen by 
considering the behavior of the function $(£"), which is monotonous and descending in the in- 
tervals (—1, 0) and (0, 1), with a discontinuity at ii^ = 0. In fact, $(£^) +ni7r changes sign in the 
vicinity of Ei, since by definition rii is an integer such that ^{Ei) + ni7r = 0. For +nj7r < 
and E — h > 0, as well as for + riiir > and E — h < 0, the nf^ term in the third sum 

in Eq. (|25|) is missing due to the step function 0(^4). Therefore, the sum is discontinuous at 
E = Ei, and the third sum over n in Eq. (ESl) can be written as a sum over E^ in Eq. (pISl). 



Appendix B 

In this Appendix we derive the conditions for the appearence of the spin polarized states from 
Eq. (|13D , written in the form = nvr, n = 0, ±1, ±2,..., where = arctanf2/£' + 

0/2 — d{E =1= h) for spin down (up) orientation. Let us denote the greatest number n for which 
^{E = —1) > mr by uq = ^{E = — l)/7r , where the symbol [A] represents the greatest integer 
less than A. 

Since ^{E) is monotonously decreasing in the energy intervals (—1,0) and (0,1) (with a 
discontinuity at -E = 0), there would be no solution in the negative energies interval, if 
^{E = 0~) > uqtt, or exphcitly if 



(f)/2±dh-'K 



> 



TT 



(j)/2±dh-d 



TT 



(26) 



One can check that the condition (^), which means that there is no spin down (up) bound 
states in the negative energies interval (—1,0), is satisfied provided that 



0/2 ±dhe ((2A; + l)7r/2, {k + l)n - d}j 



(27) 



Therefore, to have a spin up polarized ground state, 0/2 + rf/i must be inside and 0/2 — dh 
outside the above (open) interval, and vice versa for the spin down polarized ground state. 
Similar conditions hold in the positive energies interval (0, 1). There is no spin down (up) states 
provided that 



0/2 ±dh-d 



> 



TT 



^/2±dh + 7i/2 



TT 



f28) 



1. e. 



0/2 ±dhe (kix + J, {2k + l)7r/2) . 



(29) 
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For d > TT /2 neither of the conditions ( ]27| ) nor p9| ) is satisfied, so there is no spin polarized 
states. Note that the condition for ZES is obtained when 0/2 ± dh = {2k + l)vr/2, i. e. at the 
only common point of two intervals in (|2^) and (pUj). 
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Figure 1: 

Schematic illustration of the weak link. 



Figure 2: 

Normalized supercurrent at zero temperatures for (a) S/N/S, and (b) S/F/S 
weak links. In both cases the reduced thickness of the barrier is 2d/C,o = 1- 
For S/F/S two examples are presented: a "0-contact", (j)gs = 0, with h/A = 
1.5 (dotted curve), and a "vr-contact", (pgs = vr, with h/A = 4.7 (solid curve). 



Figure 3: 

Densities of states in the barrier of S/N/S weak link with 2d/C,o = 1- (a) 
PDOS for ip = 0, (b) DOS averaged over the spherical Fermi surface for 
= (solid curves), and for = vr (dashed curves). 



Figure 4: 

DOS in the ferromagnetic barrier of S/F/S weak link with 2d/C,o = 1 and 
h/A = 0.6: (a) ground state, = 0, (b) = vr. Spin splitting: a =t (solid 
curves) and a =| (dotted curves). 



Figure 5: 

DOS in the ferromagnetic barrier of S/F/S weak link with 2d/C,o = 1 aiid 
h/A = n^/2: (a) = 0, (b) ground state, = vr. Spin splitting: a =| (solid 
curves) and a =| (dotted curves). Note that for = ZES appear for both 
spin orientations in accordance with Eq. ([16D 



Figure 6: 

DOS in the ferromagnetic barrier of S/F/S weak link with 2d/^Q = 1 and 
h/A = 1.5: (a) ground state, = 0, (b) = tt. Spin splitting: a =] (solid 
curves) and a =| (dotted curves). 
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